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NOMENCLATURE 


B  -  General  linkage  pin  joint  or  points  in  links  m,n 

m,n 

C„        -  Viscous  damping  coefficient  of  dashpot  attached  to  link  Z 

Z 

and  ground  , 

* 
C  -  Equivalent  viscous  damping  coefficient  associated  with  the 

ij 

i    input  link  due  to  a  unit  angular  velocity  of  j    input 
E  -  Coordinate  point  denoting  center  of  gravity  of  system  link 

F  -  General  output  coordinate  point 


F  -  External  force  acting  through  general  system  point  E 

e 

* 
r 
e/i 


*  th 

F  -  Equivalent  force  acting  on  i    system  input  due  to  a  unit 


external  force  at  point  E 

.g         -  Velocity  influence  coefficient  of  link  I   with  respect  to 
X  1 

input  link  i 

-G.        -  Row  vector  of  velocity  influence  coefficients 

X  i 

.h         -  Acceleration  influence  coefficient  of  link  Z   with  respect 
Z   ij 

th      th 
to  the  i    and  j    input  links 

„H  .       -  Square  matrix  of  acceleration  influence  coefficients 
I   ij 

i  -  Denotes  input  link  or  generalized  coordinate 

I  .         -  Effective  moment  of  inertia  of  link  Z   taken  about  its 

center  of  gravity 

* 
I.  .        -  Equivalent  moment  of  inertia  term 


viii 


j  -  Denotes  input  link  counter 

k  -  Corresponding  time  counter 

k  -  Corresponding  position  counter 

K«  -  Effective  spring  constant  associated  with  a  spring  attached 

between  link  Z   and  ground 

K,..  -  Equivalent  spring  constant  of  K.  with  respect  to  system 

input  i 

Z  -  Denotes  general  system  link 

L  -  Lagrangian  of  the  equivalent  mass  system 

m  -  Total  number  of  general  system  outputs 

M«  -  Effective  mass  of  link  Z   at  center  of  gravity 

n  -  Denotes  total  number  of  system  inputs 

N  -  Denotes  total  number  of  system  links 

P  -  Total  system  inertia  power 

*  th 

P        -id    equivalent  inertia  power  coefficient  with  respect 

th 
to  the  r   input 

q  -  Denotes  input  link  counter 

r  -  Denotes  input  link  counter 

S  -  Summation  representing  the  time  integral  of  the  Lagrangian 

S„         -  Summation  representing  the  time  integral  of  the  virtual 

work  of  the  nonconservative  forces 

T-         -  General  system  external  torque  acting  on  link  Z 

*  th 

T,       ..  -  Equivalent  torque  acting  at  i    input  link  resulting  from 
*,e,d,s/i 

unit  torques  on  link  Z,    unit  forces  through  system  point  E, 
unit  velocities  on  equivalent  viscous  damper,  and  unit 
displacements  on  equivalent  springs 
ix 


T  -  Total  equivalent  torques  acting  on  1    system  input  link 

i 

v  -  Linear  velocity  of  general  system  point  E 

e 

w  -  Weight  of  link  I   acting  through  center  of  gravity 

x  -  Denotes  x-coordinate  of  system  input 

X  -  Denotes  x-coordinate  of  system  output 

y  -  Denotes  y-coordinate  of  system  input 

Y  -  Denotes  y-coordinate  of  system  output 

a.  -  Angular  acceleration  of  link  i 

6  -  Variation  of  some  parameter 

A  -  Finite  increment 

T)  -  Denotes  angles 

X  -   Undetermined  Lagrangian  multiplier 

cp.         -  Angle  of  i    system  input  link 

.th 
§  -  Angle  of  *   system  output  link 

0).         -  Angular  velocity  of  link  i 

J   J        -  Denotes  column  vector 

-  Denotes  differentiation  with  respect  to  time 

-  Denotes  transpose  of  matrix 

S  -  Denotes  partial  differentiation 

Sign  Convention 

Right-hand  Cartesian  coordinate  system 

Angles  measured  positive  ccw  from  positive  x-axis 
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The  holonomic  constraints  associated  with  complex,  multiple 
input  linkage  systems  complicate  the  procedures  and  methods  used  in 
determining  their  dynamic  response.   Large  systems  of  nonlinear, 
second-order  differential  equations,  requiring  additional  algebraic 
equations  of  constraint,  occur  as  a  result  of  these  constraints. 
Double  iteration  algorithms,  which  are  both  time-consuming  and  subject 
to  error,  are  necessary  to  integrate  numerically  these  differential 
equations  of  motion. 

In  this  dissertation  the  concepts  of  kinematic  influence  coef- 
ficients of  complex,  planar,  rigid  link  mechanisms  with  multiple  inputs 
are  developed  and  utilized  to  eliminate  the  holonomic  constraints  asso- 
ciated with  such  systems.   Kinematic  influence  coefficients  associated 
with  series  and  parallel  linkage  combinations  are  developed,  based  on 
the  addition  of  Assur  groups  (dyads,  tetrads  and  more  complex  groups) 
to  the  basic  system  group. 

xi 


These  complex,  multiple  input  linkage  systems  are  then  reduced 
to  coupled  equivalent  mass  systems  acted  upon  by  variable  rate  springs, 
variable  coefficient  viscous  dampers,  and  equivalent  external  forces 
and  torques.   The  holonomic  constraints  associated  with  the  original 
system  are  eliminated,  thus  leaving  the  equivalent  mass  system  free  of 
all  such  constraints.   The  number  of  generalized  coordinates  required 
to  describe  the  motion  of  the  equivalent  system  now  equals  the  number 
of  independent  system  inputs. 

The  differential  equations  of  motion  describing  the  system's 
dynamical  behavior  can  then  be  determined  by  established  methods  and 
put  in  a  suitable  form  for  numerical  integration. 


CHAPTER  I 
LITERATURE  SURVEY 

A  substantial  survey  of  the  literature  representing  the  state 
of  the  art  in  dynamic  response  analysis  of  constrained  mechanisms 
appears  in  ri].   Virtually  all  systems  treated  in  the  literature  possess 
a  single  degree  of  freedom.   Furthermore,  these  systems  are  analyzed  for 
their  output  dynamic  state.   Few  researchers  have  attempted  to  deter- 
mine the  input  dynamic  state  of  complex  mechanisms  with  multiple  inputs 
(degrees  of  freedom)  to  forces  and  torques.   Those  who  have,  invariably 
give  examples  which  are  constrained  four-link  mechanisms,  thus  leaving 
their  methods  essentially  untested.   This  emphasizes  the  need  for  a 
formalized  procedure  to  model  complex  systems  mathematically,  utilizing 
direct  methods  for  determining  their  dynamic  response  to  known  forces, 
torques  and  energy  crossing  the  system  boundaries.   This  will  eliminate 
iterative  procedures  necessary  in  solving  the  differential  equations 
of  motion. 

Dynamic  State  Analysis 

As  stated  in  [1]  "dynamic  state"  implies  that  the  velocity  and 
acceleration  of  each  point  of  every  link  of  the  system  are  piecewise 
continuous,  diff erentiable  functions  of  the  input  characteristics.   In 
1957  Modrey  [2]  developed  a  graphical  method  whereby  the  velocities 


and  accelerations  of  the  links  of  a  system  of  higher-order-complexity 
(one  which  cannot  be  analyzed  as  a  system  of  four-link  mechanisms  con- 
nected in  series)  could  be  determined,  using  velocity  and  acceleration 
influence  coefficients  obtained  through  a  "zero-relax"  procedure.   In 
a  discussion  of  this  paper,  T.  P.  Goodman  showed  the  analytical  equiv- 
alent to  Modrey's  method  and  proved  the  linearity  property  necessary 
for  superposition  with  respect  to  single  degree  of  freedom  systems. 
Other  graphical  methods  existed  for  solving  the  dynamic  state  question 
of  complex  mechanisms;  namely,  Hall    and  Alt's  method,  Carter's  method 
and  the  method  of  normal  accelerations.   All  of  these  methods  assume 
a  known  geometric  configuration.   That  is,  the  point  path  is  known,  and 
the  dynamic  state  is  determined  based  on  this  known  configuration. 

Displacement  Analysis 

Various  methods  have  been  developed  for  the  displacement  analysis 
of  mechanisms  containing  more  than  one  vector  loop.   Hain  [3]  shows  how 
to  determine  the  input-output  position  information  of  six-  and  eight-bar 
mechanisms  which  do  not  contain  a  basic  four-bar  at  ground  but  do  contain 
an  internal  four-link  loop.   This  is  accomplished  by  inverting  the  mechan- 
ism, solving  the  position  information  and  then  reinverting.   However, 
most  methods  rely  on  harmonic  analysis,  including  those  by  Meyer  zur 
Capellen  [4] ,  Denavit  and  Hasson  [5] ,  Flory  and  Wolford  [6] ,    Markus  and 
Tomas  [7],  and  the  Romanian  school  [8,9].   They  obtain  the  harmonics  by 
approximate  numerical  means  within  a  prescribed  error  criterion. 
Crossley  and  Seshachar  [10]  analyze  the  displacement  of  planar  Assur 


groups  by  an  iterative  algebraic  method  which  allows  half  the  unknowns 
from  the  matrix  of  complex  numbers  to  be  found  first  with  the  other 
half  found  later.   All  of  these  methods  rely  on  iteration  or  approx- 
imations because  of  the  complex  loop  equations  resulting  from  the 
geometry.   However,  they  yield  results  which  can  be  made  as  accurate 
as  desired.   Digital  computation  capability  increases  the  desirability 
of  these  methods. 

Dynamic  Response  Analysis 

Time  response  of  mechanical  systems,  spatial  as  well  as  planar, 
by  use  of  Lagrangian  mechanics  has  been  advanced  during  the  most  recent 
years.   Chace  [11]  uses  relative  coordinates  in  determining  the  dynamic 
response  of  multiple  degree  of  freedom  spatial  mechanisms.   He  utilizes 
Lagrangian  multipliers  to  account  for  the  physical  (geometric)  con- 
straints.  Smith  [12]  employs  this  same  technique  in  analyzing  the 
reaction  forces  in  generalized  machine  systems.   In  both  [11]  and  [12] 
the  examples  are  planar  single  degree  of  freedom  four-link  mechanisms, 
leaving  the  general  nature  of  the  method  untested. 

Uicker  [13]  and  Carson  and  Trummel  [14]  apply  matrix  notation 
and  methods  to  the  kinematics  problem.   By  applying  Lagrangian  techniques, 
they  develop  the  system  differential  equations  based  on  the  vector  loop 
displacement  equations.   This  method  has  been  the  most  useful  and  power- 
ful for  analyzing  the  dynamic  response  of  complex  multiple  degree  of 
freedom  systems  to  date.   As  with  other  methods,  it  has  its  weaknesses. 
A  dual  iteration  algorithm  results,  which  is  ideal  for  digital  computa- 
tion but  is  time-consuming. 


The  works  of  Wittenbauer  [15],  Federhofer  [16],  and  Beyer  [17] 
are  significant  in  the  development  of  equivalent  mass  systems  and 
appear  to  be  the  most  direct  antecedents  of  the  concepts  developed  in 
this  dissertation.   Wittenbauer  and  Federhofer  introduce  the  concepts 
of  reduced  mass  in  terms  of  general  mass  content  and  velocity  ratios. 
Beyer  effectively  summarized  this  work  for  the  modern  reader  in  terms 
of  the  time  response  problem.   The  authors  in  references  [11-14]  ignore 
this  property,  preferring  to  treat  these  systems  as  ones  with  large 
numbers  of  generalized  coordinates,  coupled  by  algebraic  equations  of 
constraint.   This  property  necessitates  the  use  of  the  Lagrangian 
X-method  in  treating  the  holonomic  auxiliary  conditions,  replacing  the 
kinematical  constraints  by  the  forces  necessary  in  maintaining  them. 
The  \' s,  then,  indicate  the  degree  to  which  the  constraints  are  violated. 

The  need  for  a  more  concise  system  formulation,  allowing  one  to 
determine  the  dynamic  response  in  a  more  precise  and  direct  manner,  is 
paramount.   This  will  be  accomplished  in  the  following  chapters  by 
developing  the  concepts  of  kinematic  influence  coefficients  of  complex, 
planar,  rigid-link  mechanisms  with  multiple  inputs  (degrees  of  freedom). 
These,  in  effect,  account  for  the  holonomic  constraints  on  the  system, 
resulting  in  coupled,  equivalent  mass  systems  requiring  only  as  many 
generalized  coordinates  as  its  degrees  of  freedom.   These  equivalent 
systems  can  then  be  uncoupled  by  standard  methods,  allowing  their  time 
response  to  be  determined  by  direct  methods  such  as  that  developed  in 
[18]. 


CHAPTER  II 


KINEMATIC  INFLUENCE  COEFFICIENTS  OF  COMPLEX  MECHANISMS 
CONSTRUCTED  FROM  ASSUR  GROUPS 


Pelecudi  [19,20]  defines  kinematic  multipoles  as  kinematic 
chains  interpreted  as  lines  for  transmitting  motion  information  (posi- 
tion, velocity,  acceleration)  of  rigid  links  from  an  input  to  an  output 
or,  more  generally,  from  one  point  in  the  chain  to  another.   This  is 
accomplished  by  motion  transforms  called  kinematic  ratios,  velocity  and 
acceleration  ratios,  or,  more  precisely,  kinematic  influence  coefficients. 
These  are  well  established  for  the  case  of  four-link  mechanisms  and  cams 
in  [1,21-25].   Kinematic  chains  receive  information  from  input  poles  and 
transmit  it  to  output  poles.   The  generality  of  this  concept  enables 
one  to  construct  highly  complex  mechanisms  consisting  of  groups  of 
links  known  as  Assur  groups.   The  degree  of  mobility  of  the  final  sys- 
tem can  therefore  be  controlled  or  determined  by  the  mobility  associated 
with  each  group  added  to  the  basic  kinematic  chain.   Two  types  of  chains 
are  most  important  in  the  construction  of  complex  mechanisms. 

(1)  Structural  groups  having  a  degree  of  mobility  (freedom) 
f  ^  0,  forming  pa^dLve.  multipoles.   A  group  may  be 
overconstrained  (f  <  0)  to  form  a  mechanism  with  a 
degree  of  mobility  which  is  less  than  the  basic  chain. 


(2)  Mechanisms  with  a  fixed  link  having  a  degree  of  freedom 
f  £  1,  forming  active,    multipoles. 

Thus,  one  is  free  to  construct  mechanisms  from  basic  chains 
whose  degree  of  freedom  remains  invariant  under  the  addition  of  aotcve. 
and  paAdJ-Ve.   multipoles.   This  points  out  the  importance  of  the  struc- 
tural group  in  the  construction  of  mechanisms  of  high  complexity. 

For  the  purpose  of  analysis  it  is  desirable  to  reduce  these 
systems  to  their  basic  kinematic  groups  in  order  to  derive  the  input- 
output  transformation  of  position,  velocity,  and  acceleration.   Systems 
possessing  these  transformations  can  be  reduced  to  coupled  equivalent 
mass  systems.   This  coupling  will  be  quadratic  and  will  depend  on 
velocity  influence  coefficients. 

It  is  necessary  to  discuss  the  types  of  kinematic  groups  avail- 
able for  constructing  complex  mechanisms  of  desired  mobility. 

Assur  groups  [26]  are  defined  in  terms  of  kinematic  groups 
which  form  structures  when  the  outer  joints  are  fixed  after  being 
separated  from  the  input  links.   The  binary  group  is  the  simplest  of 
the  Assur  groups.   Figure  2-1  illustrates  its  use  in  constructing  the 
four-,  five-,  and  six-bar  mechanisms.   Each  binary  group  forms  a  struc- 
ture when  it  is  separated  from  its  primary  or  intermediate  input  and 
has  its  outer  joints  fixed.   The  dyad  in  the  four-  and  five-bar  mechan- 
isms and  the  dyad  attached  to  the  intermediate  input  of  the  six-bar 
mechanism  form  dipoles  while  the  dyad  in  Figure  2-l(c)  attached  to  the 
primary  input  forms  a  tripole,  capable  of  transmitting  information  from 
the  system  input  through  the  output  pole  to  the  dipole. 


B 


jg 


(a) 


(b) 


(c) 

Figure  2-1.   Binary  Groups 


The  next  Assur  group  of  higher  complexity  is  the  tetrad   (see 
Figure  2-2 (a)).   When  the  outer  joints  A,  0  ,  and  E  are  fixed,  the 
group  becomes  a  structure  with  f  =  0.   This  group  requires  two  vector 
loop  equations  to  determine  the  angular  relationships  between  the  links. 
More  complex  Assur  groups  are  illustrated  in  Figures  2-2(b)  and  (c) . 
These  and  other  groups  are  classified  and  analyzed  for  their  displace- 
ments by  Crossley  and  Seshachar  [10]. 

Markus  and  Tomas  [7]  define  a  system  group  as  a  group  of  links 
which  has  zero  degrees  of  freedom  when  the  input  links  are  fixed.   The 
system  groups  contain  Assur  groups  but  are  more  general,  enabling  the 
authors  to  solve  for  displacement  information  of  any  system  group  by 
a  single  harmonic  analysis  algorithm.   Curtis  and  Tomas  [27]  point  out 
the  importance  of  the  binary  group  in  determining  position  information 
of  mechanisms  of  varying  complexity.   They  utilize  the  law  of  cosines 
in  developing  a  closed  form  algorithm  to  obtain  the  position  information 
of  a  large  class  of  mechanisms  consisting  of  various  binary  group  con- 
figurations. 

Kinematic  Influence  Coefficients 
of  Primary  Input  Poles 

Velocity  Influence  Coefficients 

An  implicit  functional  relationship,  f(x,y)  =  0,  describing  the 
path  of  some  general  system  point  E  (see  Figure  2-3(a))  associated  with 
the  outer  joint  of  a  binary  group,  can  be  written  based  on  known  posi- 
tion information.   In  general,  f(x,y),  g  ,  h  ,  m  ,  cp  „ ,  and  cp  ,  „  will  be 

e        e        el  jUI 

known   as   a  function  of   some   system  input  parameter   such   as  cp . .      It   is 


7777  7777 

(0 

Figure  2-2.   Assur  Groups 
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Figure  2-3.   General  System  Point 
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noted  that  g   and  h   are  defined  in  [1]  as 
e      e 


dS 
-  _   e 


'e   dcp.  ' 


d2S 
r e 


and 


m     =   slope  of  path   tangent   of   point  s   . 
e  e 

In  order   to  proceed   through   the   chain   and   determine   the   influence 

coefficients    associated  with      F(X,Y)    the  values   for  g    ,    g    ,    h    ,    h 

must  be  determined.       It    is  noted   that    (see  Figure  2-4(a)) 

dS 

v    =  g    cp    = — —  cp     , 
x       Bx  Yi        dcp.    Yi    ' 

(2-1) 

dS 

ey  • 

Vy  "   gy  9i  "   dcp.   9i    ' 
and 


d     r^,         si         df  df 

dt  <*<*.*>}  =  d^Vx+  d?V3 


This   gives 

i.    . 

df 

"     6f 

df   -  df   -  .„   „ 

=   -r-g      +X"g      =0-  (2-2) 

OX      x         Oy      y 


(2-3) 

37 
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The  expressions  for  m  and  g  are  given  by 

m  =  — -£-  ,  (2-4) 

e    — 


and 


=y? 


+  g2    .  (2-5) 

y 


Since  g  is  known,  then 
e 


ge 
gv  =  ,  (2-6) 


and 


ST 


l     =   m  g  (2-7) 

y    e   x 


Acceleration  Influence  Coefficients 

Now  consider  h  and  h   in  terms  of  m  ,  g  and  h  . 
x      y  e   e      e 

It  is  noted  that 

h   =  a*  , 
x    x 

—     t 

h  =  a   ,  for  cp.  =  1.0,   Cp.  =  0. 
y    y       i        i 

Define  7]  as  follows  (see  Figure  2-3(b)) 

£ 

i 


(2-8) 


T]  =  tan"1  {  -£  }  .  (2-9) 


For  f(x,y)  =  0,  a  closed  path  of  point  E  in  a  general  system  link, 

g  is  always  positive.   Then  h   is  positive  if  directed  in  the  same 
e  e 

sense  as  g  ,  negative  if  directed  in  the  opposite  direction. 
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Therefore, 

h  =  ii  cos  (7])  ,  ^ 
X    6  (2-10) 

h  =  h   sin  (T))  . 
y    e 

Now  consider  a  dyad  as  a  dipole  or  binary  group  in  the  Assurian 
sense  with  points  E   and  E  having  general  point  paths  defined  by 
implicit  functions   f  (x  ,y  )  =  0  and  f_(x  ,y  )  =  0  (see  Figure  2-4(a)). 
Then, 


Cl    /l  +  m-i 


% 


*>    J 


S„     =  m   g^ 


(2-11) 


1  +  m^         y2    e2   X2 

e„ 


and 


h   =  h   cos  (7]  )  ,   h   =  h   sin  (T]  )  , 
x,     e,        1      y.     e..        1 


h   =  h   cos  (71)  ,   h   =  h   sin  (7)  ) 
X2    e2       2      y2    62       * 


(2-12) 


This  can  be  generalized  immediately  to  all  classes  of  Assur  groups 
(see  Figures  2-4(a)  and  (b))  as 
g„ 


i   i 


*'    / 


i T~  ~y.  e.   x. 

1  +  mz  i 


i 

(2-13) 


h   =  h   cos  (71)  ,   h   =  h   sin  (7]) 
x     e,       i      ye.       i 
i     i  ii 
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(x,,y,) 


(x2,y2) 


(a) 


(x-.y.) 


i"i 


(b) 

Figure  2-4.   Point  Paths  of  Assur  Groups 
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Kinematic  Influence  Coefficients 
of  Intermediate  Input  Poles 

Let  E   and  E  be  the  inputs  and  F  be  the  output  of  the  binary 
group  in  Figure  2-4(a).   Point  F  is  an  internal  joint  to  links  i   and 
l+l.     However,   the  following  derivation  holdi  for  output  poles  as 
general  points  in  either  of  the  links  i   or  X+l.   The  point  path  of 
F  can  be  expressed  as  an  implicit  function 

f(X,Y)  =  0,  (2-14) 

where  X  and  Y  are  functions  of  x  ,  y  ,  x  ,  and  y  .   Relation  (2-14)  may 

1    1    A  A 

be  expressed  by  the  projections  of  the  vector  loop  equation  for  the 
dyad, 


E.F  +  FE   +  E  E  =  0,  (2-15) 

1  £t  £      \ 


fx(WX'yi'y2'Y)  =  °' 
and  (2-16) 

fy(Xl'VX'yi'y2'Y)  =  °- 


Velocity  Influence  Coefficients 

Since  X,Y  are  functions  of  x   y  ,  x   and  y  ,  then  they  are 
differentiable  with  respect  to  them.   The  velocity  of  point  F  is  given 


VF  =  vx  +  vy  (2-17) 
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where 


**  ■  £  { 


dx 


-v     J    , 


.   dx        x .         oy .       y 
i=l  i         1  'i      Ji 


dX 

d7~ 


(2-18) 


Ef    dY  Sy 

t^vx.+  3y~V 
1=1  1        1  11 


The  partial  derivatives  in  Equations  (2-18)  represent  the  velocity 

influence  coefficients  for  multiple  input  linkage  systems.   They  are 

defined  as  follows  in  order  to  remain  consistent  with  previous  literature. 

Let 

A  dx 
^x.  ~  ox"  ' 


yV 


A  dY 


(2-19) 


Substituting  these  into  Equations  (2-18)  gives 
VX=  t    {x^x.  Vx.  +  Xiy   VyJ.  ' 


i=l 
2 


Vy  =   E   {A.  Vx.  +  yiy.  Vy  } 

1=1     ii     'i  Ji 


(2-20) 


Equations  (2-20)  can  be  put  into  compound  matrix  form  as 


XGx. 


YGx. 


xGy, 


yV 


(2-21) 


Similarly  for  angular  output  information 
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M" 


A.     1      *V 

II  J  1 


(2-22) 


where 


XGx. 


"  Bx  a  x" 

dx  dx„ 


A. 


"as,        sa; 
3xi        ^ 


(2-23) 


Acceleration  Influence  Coefficients 


Differentiating  Equations  (2-18)  with  respect  to  time  yields 


2        2 


=  y    y  {  sx_  v     v     +  2    ? 

.*-*      h*.     Ldx .  ox .      x .      x .  ox 


d'x 


"X         (->      Cu     idx  dx        x.      x.  dx.oyT     x.      y. 

i=l  j=l  ijij  iJiJ 


z 

1         y      f  SX       t  dX       t   \ 

V  J          A*     Idx        x  dy        y.J 

i   *j      'i      'j           1=1          li  'i        i 


d'x 

+    5 S v 


(2-24) 


and 


2        2 


s  y    y  {,3,Y   v    v    +  2    « 

*■*      .*-*,     ldx .  dx .      x .      x .  dx 


a2v 


— T V  V 

dy        x .      y  . 

i=l  j  =  l  ijij  i      J         i        J 


^2 
d  Y 


v        v      h  +     Y*     "W —  a        +   -r —  a      r    .  (2-25) 

y.      y.J  ,     Wx.      x.         dy.      y.J 


f   dy    dy        y.      y.J        /-»     ldx.      x.         dy .      y. 
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Equations  (2-24)  and  (2-25)  can  be  written  in  compound  matrix  form  and 
generalized  to  n  inputs  as 


aX  "  'Vx' 


i 

|v    |    +   2|v    | 

XHx.y. 

|v    1 

X.  X  . 

1      X                   '      X1 

y 

1  J- 

i  J_ 

"l  * 

_                __ 

a 

1      X 

f 

|v     |     + 

VG 

Jg 

y-y  • 
1  J_ 

1        yl 

X   x. 

l 

|       XyiJ 

1    t 

• 

L  y 

(2-26) 


YHx.x. 
i  J 


fv  |  +  2|v  T  T^    "||v| 
1  x'     '  x1    Y  x.y .   ■  y1 

L    i   JJ 


+  Iv  I     H      [v  I  +   G 
1  y'    Y  y±yj   '  y1    Y  x. 


yV 


(2-27) 


where 


i,j  =  1,2,... ,n, 


XHx.x. 


xVx.    X11. 


"  X1X2     *  XlXn 


„h       __h     •••   h 

A    X2X1  A    X2X2  X    X2Xn 


„h  h  • •  •      h 

Xx  x„         X   x  x„  Xxx 

n   1  n  2  n  n 


(2-28) 


and 


xVx     5 
i   J 

32x 

dx   dx 

< 

1 

K 

= 

Xgx,    '     '    Xgx     XgyH    '     '    Xgy 
1  n        Jl  Ji 


(2-29) 


(2-30) 
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Corresponding  expressions  relating  angular  input  to  angular  or 
linear  output  can  be  derived.   Since  most  kinematic  chains  consist  of 
lower  pair  connections,  the  remaining  derivations  will  be  based  on 
angular  inputs  to  both  linear  and  angular  outputs.   Complex  kinematic 
ehaini  constructed  from  basic  chains  with  invariant  mobility  arc  dcvclopod 
by  adding  Assur  groups  Lo  output,  polos  of  the  established  chain  of  pre- 
determined mobility.   If  the  mobility  does  not  remain  invariant,  then 
the  basic  kinematic  chain  with  its  known  input-output  position  relation- 
ships is  destroyed  and  the  point  paths  are  no  longer  predetermined  with 
respect  to  the  original  system  group  inputs.   A  new  system  group  must  be 
established  and  new  input-output  relationships  developed. 

Total  Input-Output  Kinematic 
Influence  Coefficients 

If  the  inputs  to  a  linkage  system  group  are  angular,  then  the 

intermediate  linear  inputs  (x. ,y.)  and  their  dynamic  states  (vx. ,Vy.) 

t   t 
and  (ax. , ay . )  can  be  expressed  as  functions  of  the  input  position  param- 
eters CD  's  and  their  dynamic  states  uu  's  and  a   's.   The  derivation  of 

q  q      q 

the  total  input-output  transformation  is  now  carried  out. 

Velocity  Influence  Coefficients 

Consider  the  system  group  in  Figure  2-5.   The  expression  for 

the  velocity  components  of  point  F  with  respect  to  points  E   and  E   are 

1      £ 


XGx. 


YGx. 


xGy 


vGy 


1        '   X1 


(2-31) 
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E2(x2,y2) 


Figure  2-5.   Seven-Link  System  Group 
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The  intermediate  input  data,  in  terms  of  the  primary  input  information, 
are  given  by 


x.C9 

i   q 


y.G9 

i   q 


CD   ,   i,q  -  1,2; 


where 


G 

X.  CD 

i  Yq 


8*1 

*1 

ax2 

5?I 

3X2 

and 


(2-32) 


Substituting  the  right-hand  side  of  Equation  (2-32)  into  Equation  (2-31) 
yields 


xc*. 


Y°x. 


xGy, 


vV 


G 

X  cp 

i  q 


G 


which  reduces  to 


(2-33) 


X  cp 


U)  ,     q  =  1,2. 


(2-34) 


YGcp 
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Note  that  when  q=l,  Equation  (2-34)  becomes 


'1   1 


(2-35) 


VY  =  Ygl  Wl  ' 


which  are  recognized  as  the  velocity  influence  coefficients  developed 
in  reference  [1]  for  single  input  systems. 

.Acceleration   Influence  Coefficients 


The  component  accelerations  of  point  F  (Figure  2-5)  in  terms 
of  the  intermediate  inputs  are 


LI, 


i  j 


Y  x.  x . 
i   J 


X»x 


iyj 


A 


iyj 


V  + 

y1 


L 


1  !  Xy  .  y  . 

I,  i    J 


;    v 


JL*Vj 


XGx. 


Y°x. 


xV 


yV 


1    x' 

|a*| 

(2-36) 


The  first  and  third  terms  of  Equation  (2-36)  are  recognized  as  real 
quadratic  forms  while  the  second  is  real  bilinear. 
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Substitution   of   Equations    (2-26) ,     (2-27)  ,    and    (2-32)    into 


Equation    (2-36)    gives 
(2.Nxl)  (2N   X  2n) 


(2n    x  2n) 


(2n  xn)  (n  xn)   (n  xD 


t 

av 

|ffi| 

0 

X 

t 

/ 
|     | 

aY 

0 

|cu| 

b\~] 


IISM 


J^X    X. 

1  J 


YV: 


i   J 


x.  cp 
J      q 


+  2 


'w 

0 

, 

0 

M 

"w 

0 

"     T — 

/ 

0 

w 

(2N   X  2n) 


X°x. 


YGx. 


xV 


yV 


(2n 

X2nn) 

M 

0 

0 

0 

0 

l-l 

/  J 

0 

0 

1 

jl_  _ 

0 

1 

0 

!   M 
i 

I 

0 

1 

T 

/ 

0 

0 

0 

l-l 

(2nn  xn)    (nxD 


H 

x.   9  cp 
l      q   r 


y.   V  9 


(2nxn)    (nxl)' 

|or|L         N,n  =    2. 


G 

x     Cp 

i    q 


G 

y.   9 
i    q 


(2-37) 
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The  complexity  of  Equation  (2-37)  is  reduced  by  collecting  like  terms 
and  writing  it  in  a  more  concise  and  familiar  form  as 


t"I 


XHcp  9 
q  r 


YH9  9 
q  l 


xG9 


A 


\a\ ,      q,r  =  1,2 


(2-38) 


where   ||  U)|J   is  a  compound  diagonal  matrix  given  by 


(2-39) 


and   j co  I   is  the  transpose  of  the  column  vector  |u)|. 

The  expressions  (2-34)  and  (2-38)  for  the  velocity  and  acceler- 
ation can  be  generalized  to  treat  systems  with  n  inputs  and  m  outputs. 
These  expressions  take  the  general  form 


1  Tq 


x  9 
m  q 


Yl\ 


Y  9 
m  q 


hi  f   Q  =  1,2,.  .  .  ,n, 


(2-40) 


25 


t 
Y 


„  H 

Xi  ^ 
1   q  r 


m  q  r 


1  q  r 


Y  9  9 
m  qTr 


U)  + 


*l\ 


XG9 
m  q 


h\ 


G 

Y  9 
m  q 


!«  , 


(2-41) 


q,r  =  1, . . . ,n. 

It  is  clear  that  expressions  similar  to  Equations  (2-40)  and 
(2-41)  can  be  written  for  the  angular  properties  of  the  system  links. 
For  the  purpose  of  illustration,  both  forms  will  be  used  where  appro- 
priate throughout  the  remainder  of  this  dissertation. 

The  discrete  form  [28]  of  the  kinematic  influence  coefficients 
of  velocity  and  acceleration  are  expressed  as 

,8i, 
Xgik  = 


rik 


i  k 

6  _       5 

li,k+l        I  i,k-l 

2A9. 


(2-42) 


d2$ 


*hij,k  =  br4:\  • 


i,k+l 


i,k+l 


i,k-l 


i,k-l 


/ij.k 


S  _      6  _       v  +       i 

4j,k+l        Xj,k-1        X,j,k+1        Xj,k-1 

4(A9-)2 


).  (2-43) 
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where  the  input  type  (i.e.,  angular  or  linear)  is  considered  understood 
and  its  designation  is  dropped  for  simplification  of  notation  as  follows; 


(2-44) 

Table  2-1  is  a  compact  collection  of  the  influence  coefficients 
derived  thus  far,  expressed  in  discrete  form.   Each  column  is  headed  by 
a  characteristic  set  of  input  conditions  and  each  row  is  designated  by 
a  particular  system  model.   Rows  and  columns  intersect  in  blocks  con- 
taining the  influence  coefficients  for  which  these  properties  hold. 
Blocks  1.1  -  1.3  and  3.1  -  3.3  have  been  derived  and  are  expressed  in 
discrete  form. 

The  coefficients  in  blocks  2.1  -  2.3  pertain  to  the  relative 

angular  motion  between  any  two  links  m,n  in  the  system  (where  m,n  are 

not  to  be  confused  with  m,n  denoting  the  number  of  inputs  and  outputs 

as  in  Equation  (2-41)).   The  links  are  shown  with  coincident  points  at 

B  for  purposes  of  simplicity  of  notation.   The  influence  coefficients 

in  blocks  4.1  -  4.3  apply  to  the  relative  motion  between  any  two  points 

B  ,B   in  different  links  m.n  in  the  system.   The  system  model  has  been 
m  n  •  j  j 

represented  as  a  sliding  pair  connecting  these  two  points.   It  may  be 
helpful  to  consider  this  pair  as  massless  when  it  is  not  required  for 
constraint  in  the  system.   An  example  of  when  it  would  be  required  for 
constraint  is  that  shown  in  Figure  2-6. 

The  influence  coefficients  in  Table  2-1  are  a  precise  defin- 
ition of  the  system  geometry  up  to  the  second  order.   As  such  they 
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28 


Figure  2-6.   Sliding  Pair  Constraint 
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represent  a  powerful  definition  of  the  meaning  of  multiple  input  link- 
age systems. 

These  influence  coefficients  enable  one  to  determine  the 
dynamic  state  of  every  link  in  a  general  linkage  system,  given  the 
dynamic  state  of  the  system  inputs.   They  will  be  used  in  the  follow- 
ing chapter  to  reduce  complex,  multiple  input  linkage  systems  acted 
upon  by  external  forces  and  torques,  internal  springs,  and  viscous 
dampers,  to  equivalent  mass  systems  acted  upon  by  equivalent  external 
torques,  variable  rate  springs  and  variable  coefficient  viscous  dampers. 

The  remaining  coefficients  will  be  treated  on  a  discerte  time 
basis,  since  some  of  them  will  depend  on  their  past  history.   This 
means  that  expressions  such  as  that  for  S  imply  that 

S  =  f(t)  .  (2-45) 


CHAPTER  III 


EQUIVALENT  SYSTEM  FORMULATION  IN  TERMS 
OF  KINEMATIC  INFLUENCE  COEFFICIENTS 


The  coefficients  developed  in  Chapter  II  will  now  be  utilized 
to  eliminate  the  holonomic  constraints  associated  with  complex  linkage 
systems,  reducing  them  to  coupled,  equivalent  mass  systems.   These 
same  coefficients  will  be  used  to  reduce  the  generalized  internal  and 
external  force  generators  acting  on  the  system  to  equivalent  general- 
ized torques  acting  at  the  system  inputs. 

The  Equivalent  System  Torques 

The  force  related  influence  coefficients  are  derived  by  replac- 
ing the  effect  of  external  forces  on  the  system  (i.e.,  T .,    F  )  and  inter- 
nal forces  generated  by  system  elements,  such  as  springs  and  viscous 
dampers,  with  equivalent  torques  acting  on  the  n  system  inputs. 

External  Forces  and  Torques 

If  a  linkage  system  (see  Figure  3-l(a))  is  given  a  virtual 

displacement  by  each  separate  input  from  some  given  system  position  k, 

th  . 
then  the  virtual  work  done  by  the  equivalent  torque  at  the  1    input 

must  be  equal  to  the  virtual  work  done  by  the  system  external  torque,  or 


TVi  *Jk  ■  lTi  A§Jk'  .  (3-1} 

30 
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(a) 


%%*' 


K 


Vi 


(b) 


* 


(c) 


Figure  3-1.   Equivalent  System  Elements 
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solving  for  T^/iR  and  taking  the  limit  as  ^?ik"*0  gives 


J 


'4k 


Similarly,  for  an  external  system  force  F   , 


TX/ik   Ugikj  Fek 


(3-2) 


(3-3) 


The  total  equivalent  torque  seen  at  each  input  link  i,  resulting  from 
all  system  torques  T .    ,    is  given  by 


,Si 


1&1   2gl  '  '  ^1 
1S2    2g2  '  '  1^2 


lgn   2gn  '  '  ^n 
i_   _tk  u 


Tl 

lGi 

Tl 

*M 

NGi 

\ 

Ik 

k 

k 

L-    -1 

(3-4) 


i  =  1,2 n, 

N  =  number  of  system  links. 

Equations  (3-2)  and  (3-3)  show  that  „g.,  and  g   are  the  influence  of 

Hi   ik     e  lk 

th 
a  unit  torque  or  force  applied  to  system  link  *  or  point  E  on  the  i 

system  input  link.   The  «S-.     is  the  mechanical  advantage  of  link  I   with 

th 
respect  to  the  i   input  link  and  vice  versa. 


Internal  Springs 

Perhaps  the  most  difficult  system  parameter  to  develop  in  terms 
of  an  equivalent  coefficient  is  the  equivalent  spring  constant  K«/ik- 
Suppose  a  torque  T ,      is  generated  by  a  linear  spring  K.  between  link  I 
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and  the  ground  link.   The  question  arises:   What  is  the  equivalent 
spring  constant  K,  .   of  the  system  spring  K  .?      By  giving  the  system 
a  virtual  displacement,  the  potential  energy  change  in  the  system 

spring  can  be  equated  with  the  potential  energy  in  the  equivalent 

th 
spring  associated  with  the  i   input  link  (see  Figure  3-2(b)).   This 

yields  the  following  relation 

I  KI/ik  GO  -  §  ej  KJ  •  (3-5) 

Solving  for  K,.   and  taking  the  limit  as  #p   -•  0  yields 

*vik-rj  .  lim~  (*r).  =  At*i  ■  (3_6) 


^ik 

It  is  noted  that  §»  cannot  continuously  increase  with  increas- 
ing cp.  ,  otherwise  the  spring  would  be  destroyed.   Rather,  §.  is  cyclic 
i  * 

or  quasi-cyclic  with  increasing  9..   Therefore,  for  part  of  the  cycle 

A§.  (i.e. ,  .g.)  must  be  negative.   This  implies  the  need  for  the  signum 

function  a    =  sgn  („g.)  which  changes  sign  whenever  the  flg.'s  go 
fg.  Ill  *   1 

through  zero.   The  torque  T«.   is  then  determined  from  its  past  history 
in  the  form 

T*/-,     =Ta/-1     ,   +  K,,,.     tip.,  ,  (3-7) 

£/ik  £/i,k-l  Vik        lk 


where 


•U  ■  Vu r*  (^4^)  • 
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The  equivalent  torques  at  each  system  input  resulting  from  all  system 
springs  are  given  in  matrix  form  as 


s/1 


s/2 


s/n 


4   L 


s/1 


s/2 


*  * 
s/n 


11 


k-l 


K£/l  *i 


KX/2  ^2 


Kly   A? 

*/n   n 


(3-9) 


An  alternative  approach  to  determine  the  equivalent  torque 
T..   generated  by  a  system  spring  between  link  JL   and  ground  is  to 
treat  the  torque  acting  on  link  i,   due  to  spring  K.   as  an  external 
torque  as  in  Equation  (3-2).   The  external  torque  is  given  by 

* 


TXk  =  Ki 


(§Xf  "  $Xk)  =  KX 


A§ 


if        ^£kj       "X  "  Xk  ' 
where 

$,_  =  free  length  of  spring  measured  from  the  reference  axis, 
A$.  =  deformation  of  system  spring  K-  from  free  length. 
Substituting  Equation  (3-10)  into  (3-2)  gives 


(3-10) 


T£/ik  "  KJL 


CegiJ 


A$*  . 
lk 


(3-11) 


The  equivalent  torques  at  each  system  input  resulting  from  all 
system  springs  K^  are  now  given  in  matrix  form  as 


s/1 


s/n 


-•k   L 


lGi 


N°i 


k  L 


Kl  Ml 


K   A$ 
N    N 


,  i  =  1,2,. . . ,n,  (3-12) 

N  =  number  of  system  links. 
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Inter, Ml  Viscous  Dampers 

Let  the  torque  T.     be  generated  by  a  system  dashpot  C.  between 
link  I   and  ground  (see  Figure  3-l(c)).   The  torque  generated  by  the 
dashpot  is  proportional  to  the  angular  velocity  of  the  link  attached  to 
the  dashpot ,  or 


j=l  -  J  J'k 


(3-13) 


—  th 

Using  Equation  (3-2) ,  T.      is  transferred  to  the  i   system  input  as 


T4/ik  =  UsilJ  ~^l   £    \#i  wj  )     • 


(3-14) 


th 
The  total  torque  acting  on  the  i    input  link  resulting  from  all  system 

dashpots  is 


(    * 


d/ik 


3=1 


J  J 


(3-15) 


where 


<j,k=  £(^xgi  ^)k 


(3-16) 


The  equivalent  torques  at  all  system  inputs  resulting  from  all  system 
dashpots  are  given  in  matrix  form  as 


d/1 


d/2 


d/n 


C      C 
11     12 


C21     C22 


c  -,     c  o 
nl     n2 


In 


"2n 


(3-17) 
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Table  3-1  is  a  concise  collection  of  the  preceding  equivalent 

torque  concepts.   Row  2  corresponds  directly  with  row  1.   Those  in 

row  2  are  more  general ,  since  the  generated  torques  occur  between  two 

moving  links  m  and  n  in  the  system.   The  relative  influence  coefficients 

r    correspond  to  those  in  Table  2-1,  i.e. , 
b  ik 

r.,  =   g   -  g    .  (3-18) 

b  ik   n  ik   m  ik 

The  equivalent  torques  in  row  3  refer  to  a  force  F   applied 

to  a  particular  point  or  points  in  the  system.   Since  the  direction 

angle  T)   of  this  force  relative  to  the  path  tangent  is  known  as  a 

function  of  the  input  angles  cp .  ,  then  the  appropriate  equivalent  torque 

can  be  determined.   The  equivalent  torque,  T  ...  ,  is  obtained  by  first 

e/ik 

taking  its  component  F   along  the  path  tangent  and  transferring  this 

force  to  the  i   input  link  by  the  influence  coefficient  S-k-   Here, 

the  product  g   cos  (71   )  is  the  effective  influence  of  a  unit  force 
e  ik       ek 

applied  in  the  direction  F    at  point  E  on  the  i    input  link.   As  such, 

it  corresponds  directly  with  .g„  in  row  1.   Blocks  3.2  and  3.3  follow 

a   ik 

directly  as  a  result  of  this  correspondence.   Here,  the  system  spring 
K  and  dashpot  C  are  anchored  at  some  fixed  point  0  . 

Finally,  the  equivalent  torques  in  row  4  apply  to  forces  in 

extensible  two-force  members  attached  to  links  m  and  n  at  points  B  and 

m 

B  .   The  relative  influence  coefficient  r   is  defined  in  Table  2-1, 
n  b  ik 

row  4.   This  implies  the  treatment  of  very  general  force  generators  in 
the  system. 
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Total  Equivalent  Torque 

th 
The  total  equivalent  torque  acting  upon  the  i    input  link 

is  given  as 

T*  =  T*  .,  +  T*    +  T*  ...  +  T*  ..   .  (3-19) 

ik    4/ik    e/ik    s/ik    d/ik 

If  one  or  more  of  the  input  elements  is  actually  a  sliding 

oair   then  s    represents  that  input's  reference  position  parameter 
^  ik 

* 
and  F   would  be  the  equivalent  force  acting  on  that  system  input, 
ik 


Equivalent  System  Inertias 

The  objective  of  this  section  is  to  transfer  the  effective 
mass  of  link  I   to  each  of  the  n  system  inputs  in  order  to  obtain  equiv- 
alent mass  systems  for  a  linkage  system  such  as  that  shown  in  Figure  3-2. 

Consider  the  kinetic  energy  of  a  complex,  multiple  input  system. 

k  - 1  £  { h  -4 +  5*  4 }  •  (3-20) 


(KE) 


where 


I  .   =  effective  moment  of  inertia  of  link  I   about  its  center 

Jw 

of  gravity  E, 

M.  =  effective  mass  of  link  i   located  at  its  center  of  gravity  E. 

Writing  Equation  (3-20)  in  terms  of  the  system  input  velocities  yields 

N        n  2        n  _2 

(KE)k  =  |  Z  {Ti  L E  <#ik)-ikJ  -  «i  LE  (eiik)-ikJ  }  ■ 
&=i       1=1  1=1 

(3-21) 
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n  =  3 

N  =  13 


Figure  3-2.   Complex  Multiple  Input  Linkage  System 
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Rewriting  Equation  (3-21)  in  a  more  compact  form  gives 
n    n    jj 


CKE)k  =  §    £        Z        t    {  frtWW    +   MX(egi)(eij)]  vl 
1=1     J  =  l      *=1 


where  it  is  noted  that 


)  =  { 


i6i,vrj 


1    Z=i,j;  i  =  j 

0    X=i,j;    i^j, 


k 
(3-22) 


(3-23) 


a  direct  consequence  of  the  meaning  of  independent  system  inputs. 
The  kinetic  energy  can  be  written  in  matrix  form  as 


(KE)k  =  \    |o,|k 


(nxn) 


(3-24) 


where 


*     =  E 


£=1  L- 


I/#i)(#j)  +  Veg 


.)(  g.)l 

1   ejJk 


(3-25) 


Matrix  Equation  (3-24)  is  recognized  as  being  real  quadratic  form 
since  the  [I. .]  matrix  is  real  symmetric.   Thus  it  is  seen  that  any 
complex  multiple  input  linkage  system  can  be  transformed  into  a  system 
of  coupled  equivalent  masses  with  n  degrees  of  freedom  whose  equiva- 
lent inertias  are  given  by  Equation  (3-25). 
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Two  Degrees  of  Freedom  Example 

For  the  purpose  of  illustration,  consider  the  differential 
gear  (Figure  3-3)  ,  where  the  translating  links  are  constrained  to  move 
in  the  horizontal  direction  and  rotation  of  link  Z     is  positive 
countorclockwise. 

The  linear  velocity  influence  coefficients  of  &     are 

-     SX3    _^2_ 
3gl  =  3^  =  di+d2   ' 


(3-26) 


-    SX3   _^1_ 
3g2  "  o^  -  di+d2 


and 


v  =   s:   v  +   k   v  (3-27) 

3    3B1   1    3&2   2  v 

The  angular  velocity  influence  coefficients  of  SL     are  given 
by  (see  Figures  3-4  and  3-5) 

4gl  "  ^l  -    ~   d1+d2  ' 

(3-28) 

S$4      _±_ 
4g2   3x7   +  d.+d 


and 


2      12 


^4  =  4gl  Vl  +  4g2  V2  •  (3_29) 


The  kinetic  energy  of  the  two  input  differential  gear  is 
1  ,-  2    -2    -   -.2    -2. 


KE 


=  2  i'Vl  +  M2V2  +  (M3  +  M4)V3  +  W  •  (3'30) 
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dx, 

1 

>^*-~- 

^&/\ 

< 

i / 

V — 

dx,  =  (d,+  d8){-d$4) 


Figure  3-4.   Angular  Relationship  to  Input  #1 


d*. 

_\d$4 

dx, 


dx2=  (d,  +  d2)d$4 


Figure  3-5.   Angular  Relationship  to  Input  #2 


Substituting  Equations  (3-27)  and  (3-29)  into  (3-30)  gives 


KE  =   \   {MlV2  +  U2v22   +  (M3  +  M4)(3g1vi  +  3g2v2)2 
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+  I4(4glVl  +  4g2V2}  )  ' 


(3-31) 


where 


Collecting  like  terms  and  rewriting  Equation  (3-31)  gives 

KE   =  \   {[\  +   T4(4gl)2   +    (M3  +    V^i/lv* 
+    [M2   +   ^4(4g2)2   +    (S3   +   M4)(3i2)2]v^ 
+    2[T4(4gl)(4g2)    +    (M3  +   M4)(3i1)(3i2)]viv2} ,  (3-32) 


1     r     *       2  *        2  * 

KE   =      2    ^11V1   +    ^2*2   +    2I12V1V2} 


(3-33) 


where 


hl  =   \  +   hW*  +    (M3   +   M4)(3gl)2' 


T12=    I4(4Bl)(4g2)    +    (M3+V     (3gl)(3g2)     ' 


(3-34) 


X22   =   M2+    F4(4B2)2+    (M3  +   V(3g2)2 


For   the   system   in  Figure   3-3,    Equation    (3-24)    becomes 


KE  =  2[V1      V2] 


11 


"21 


* 



12 

V 

1 

* 
22 

V2 

—          — 

(3-35) 


which  is  recognized  as  a  quadratic  form. 
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In  the  case  of  a  differential  gear  the  equivalent  inertias  are 
constant.   However,  for  a  general  linkage  system  (Figure  3-2),  the  I  's 
are  not  constant  but  a  function  of  the  system  input  parameters.   Note 
that  whenever  v  =0,  then 

KE  =  |  I*x  v\    ,  (3-36) 

and  whenever  v  =  0 ,  then 

KE  ■  \   i  V2  '  (3"37) 

which  are  simply  the  expressions  for  the  kinetic  energy  of  a  single 
input  equivalent  system. 

Inertia  coupling  terms,  I. .   i  ^  j,  appear  for  systems  with 
two  or  more  inputs,  indicating  a  relationship  between  the  I  .  S 
through  the  velocity  product  terms  cu.cu..   Quinn' s  energy  method  [29] 
becomes  invalid  for  these  types  of  systems,  since  the  distribution  of 
kinetic  energy  between  the  links  of  the  system  is  no  longer  invariant 
with  respect  to  the  input  velocities.   The  equivalent  mass  system 
described  above  remains  valid  provided  certain  modifications  are  made 
with  respect  to  the  solution  technique. 

Equivalent  System  Inertia  Power 

The  final  equivalent  system  parameter  to  be  considered  is  the 
equivalent  system  inertia  power.   The  coefficients  derived  in  this 

section  are  necessary  when  deriving  the  differential  equations  of 

* 
motion  of  complex  linkage  systems.   If  the  I.  .  s  were  constant,  as  in 

the  preceding  example,  these  terms  would  be  zero.   Although  they  appear 
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unwieldy  in  size  and  computation  for  systems  with  more  than  two  inputs, 
they  readily  lend  themselves  to  computer  computation. 

The  power  necessary  to  drive  a  linkage  system  against  its 
inertias  is  given  as 


where 


P  =  -£  (KE)  , 
dt 


KE  =  f  w'[i;Jw  . 


i,j  =  1,2,. . . ,n 


Differentiating  Equation  (3-24)  with  respect  to  time  gives 


„*   1  I  I   d 

P  -  2  H   dt 


ij 


|u>|  +  \a\   h*.   |u)|. 


(3-38) 


(3-24) 


(3-39) 


The  second  product  on  the  right  side  of  Equation  (3-39)  is  self- 
explanatory.   The  time  derivative  of  the  equivalent  inertia  matrix, 
however,  is  not  as  trivial.   Consider  a  typical  element  of   C I  - ^ 3 
given  by 

N 

Differentiating  Equation  (3-40)  with  respect  to  time  gives 


(3-40) 


n    N 


&)■■;:  ?.{' 


k   r=l 


W^V  +.W(4hir> 


where 


+  M, 


WW  +  W  <eV   J,  "W         (3"41) 


(Xgi)(ihjr)  =  cxp-  WW 
i    J  ^r 


(3-42) 
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Equation  (3-41)  can  be  written  as  the  product  of  a  row  and  col\ 


vector  as 


dt  /    [_ij  1     ij'nj 


(3-43) 


where 


ij  r 


is 
£=1  L 


^  l^L^W  +  WW  J 


(agJ(h.J  +  (  g.)(  h.  )  1  } 
e  i   ejr     ej   e  lr  J  J 


(3-44) 


The  subscript  notation  is  defined  by 


*  a  Bl  •  • 
•  -P   =  -IT1 


(3-45) 


Substituting  Equation  (3-43)  into  Equation  (3-39)  gives 

(lxn)  (nxnn)  (nnxn) 


K  •  §  K 


ll"r 


P 

21  r 


*  1 

P  1 

12  r  ' 


P  ' 
22  r  ' 


nl  r 


0P 
n2  r 


In  r 


„  P 
2n  r 


I 

1    * 

!   p 

nn  r 


-"k  L_ 


0 


t- 


l<  fci  K  ■ 


(3-46) 


where 


P 
ij  r 


is  a  compound  nxnn  matrix  whose  submatrices  are  lxn 


row  vectors. 
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Two  Input  Example 


Consider  the  system  in  Figure  3-6  whose  kinetic  energy  is 


given  by 


(KE)k  =  -  [o,1  u,2]k 


11 


21 


12 


22 


->k 


(3-47) 


where 


ij.k 


s 

=  E 


X=l  L. 


jfii  #j  +  M-e  £ 


i  *]k 


(3-48) 


P       P 
11  1     11  2 


P        P 
12  1     12  2 


Differentiating  Equation  (3-47)  with  respect  to  time  gives 

lo,  ,  o  " 

1  I 

w9  !  o 
z  i 

o    !  w 
i 

0   i  (0„ 


k   2  L  1  2Jk 


21P1 


P 
21  2 


P        P 
22  1     22  2 


ID 

L2Jk 


-Ik 


L  1  2Jk 


11 


'21 


* 

12 

V 

* 

10 

L  2 

22^ 

where 


8  _ 


(3-49) 


(3-50) 


Table  3-2  represents  the  concepts  of  equivalent  system 
inertias  and  equivalent  inertia  power  for  complex,  multiple  input 
linkage  systems.   With  this  formulation  the  holonomic  constraints  are 
eliminated  and  the  equations  of  motion  derived,  utilizing  only  those 
generalized  coordinates  associated  with  the  independent  inputs. 
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Figure  3-6.   Complex  Two  Input  Linkage  System 
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CHAPTER  IV 
TIME  RESPONSE  OF  EQUIVALENT  MASS  SYSTEMS 

Complex,  multiple  input  linkage  systems  present  a  unique  problem 
with  respect  to  determining  their  dynamic  response.   Methods  for  deter- 
mining the  time  response  of  single  input  systems  are  well  established, 
either  in  terms  of  Lagrangian  mechanics  [11-14]  or  equivalent  mass  sys- 
tems.  Each  method  has  its  strengths  and  weaknesses.   The  method  chosen 
to  solve  any  given  problem  should  be  based  on  need  and  solution  form. 
It  is  possible  to  obtain  phase  plane  solutions  for  single  input  systems 
by  using  techniques  in  [1]  if  the  input  velocity  is  never  less  than 
zero.   With  slight  modifications  in  the  predictor  equation,  negative 
velocities  can  also  be  treated.   Here  the  independent  variable  is  the 
input  angle,  cp . ,  rather  than  time,  a  result  of  the  cyclic  nature  of  most 
mechanisms.   However,  a  large  class  of  problems  does  not  meet  this 
requirement.   Hence,  a  more  general  method  must  be  employed. 

As  pointed  out  earlier,  the  energy  distribution  method  is 
invalid  for  multiple  input  systems.   Similarly,  the  kinetic  energy 
and  power  concepts  utilizing  the  equivalent  mass  and  force  system 
developed  in  [1]  fail  because  they  yield  only  one  second-order  differ- 
ential equation  in  n  independent  variables.   These  failings  are  a 
result  of  the  geometric  constraints  imposed  on  the  system. 
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Two  methods  will  be  discussed  in  this  chapter.   The  differ- 
ential equations  of  motion  for  a  two  input  system  will  be  derived, 
using  both  methods.   The  complete  set  of  first-order  difference  equa- 
tions will  be  derived  for  only  the  second  method  for  reasons  which 
will  be  explained  later.   The  methods  used  to  derive  the  differential 
equations  and  to  solve  them  are  not  new.   The  kinematic  influence 
coefficients  and  coupled  equivalent  systems  (inertias)  described  in 
Chapters  II  and  III  are  new  and  unique,  enabling  established  methods 
to  be  used  in  deriving  and  solving  the  differential  equations  of  motion. 

Lagrange's  Method 

A  system  of  n  second-order  differential  equations  is  derived 
from  Lagrange's  equations  as 

d  /dL  \    dL     * 
St  b-)-<5T=  V   Ci=1.2,...,n)  (4-D 

o^i     Ti 

where 

cp.  -  generalized  coordinates, 

cp.  -  generalized  velocities, 

L   -  Lagrangian  of  the  coupled  equivalent  mass  system, 

T.  -  equivalent  nonpotential  forces  and  torques. 

The  n  second-order  differential  equations  resulting  from 
Equation  (4-1)  are  nonlinear,  coupled,  nonhomogeneous ,  and  contain 
variable  coefficients.   It  is  necessary  to  reduce  this  set  of  equations 
to  2n  first-order  differential  equations  in  order  to  solve  them  on 
a  digital  computer.   Although  this  reduction  is  entirely  possible,  it 
is  often  tedious,  yielding  unwieldy  equations.   The  variable  coefficients 
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are  functions  of  geometry  and  are  hence  implicit  functions  of  time. 
This  characteristic  requires  the  determination  of  the  equivalent 
inertia  power  coefficients  of  the  system,  terms  which  are  very  diffi- 
cult to  obtain  (see  Chapter  III)  and  only  adds  to  the  complexity  of 
the  problem. 

Example 

The  five-bar  mechanism  in  Figure  4-1  is  a  representative  two 

input  system.   The  generalized  coordinates  are  cp  and  cp  •  a  helical 

spring  and  viscous  damper  are  attached  between  links  3  and  4  at  pin 

joint  B,  while  external  torques  T  and  T  are  applied  to  links  1  and 

1       £ 

2,  respectively. 

The  Lagrangian  for  this  system  is  given  as 

L  =  (KE  -  PE) ,  (4-2) 

and   the  equivalent   nonconservative  torques   are 

V"S4bri    Cbri  *1   +   br2  ^2]    +T1' 


(4-3) 


T2=-C34br2  Cbri  *1  +  br2  *2]  +T2' 


where 


2    2 
KE  =  -  T    V     I*.  cp.  cp.  ,  (4-4) 


1  —  2 

PE  =  -  K    (CD     -(D).  (4-5) 

2  34  y34f   y34 


Substituting  Equations  (4-3) ,  (4-4) ,  and  (4-5)  into 
Equation  (4-1)  and  performing  the  designated  differentiation  with 
respect  to  the  generalized  coordinates  and  velocities  yields 
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7777 


Figure   4-1.      Two    Input   System 
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dT  {   *11  *1   +    J12  »2  )   +   K34    <*34f    "  »34)(bri> 


=   T1    "C34    [(br?^l+    WWV     '  (4"6) 


3?  {   X12  *l  +    T22  ^J  +  K34    (*34f    "  CP34)(br2) 


=   T2    ~C34    [(bri)(br2^1+    V^V    •  (4"?) 


Completing   the   differentiation  with   respect   to   time  gives 

*    •  *  +••  *    •  2  *••  %•• 

11P1  *1  +    (11P2   +    12V    V2   +    12P2  *2   +    Xll  »l   +    X12  *2    "   Tl 

+   ?34(CP34f    "  CP34)(bri)    -    S4  [V^l   +    (bri)(br2)^2]   ' 


(4-8) 


21*1  *\  +    (21P2   +    22VV2   +    22P2  *2   +    J12  *1   +    J22  ^2    "  T2 

+  ?34(CP34f    "  CP34)(br2)    -    Si    W^l   +    VJ* J  ' 

(4-9) 

Equations  (4-8)  and  (4-9)  are  nonlinear  in  the  velocities  and 
are  unwieldy  in  their  present  form.   In  addition,  they  should  be  reduced 

to  four  first-order  equations  before  they  are  numerically  integrated. 

* 
The  . .P   terms  as  defined  by  Equations  (3-44)  and  (3—15)  are  not  simple 
ij  r 

expressions,  and  require  calculation  at  each  integration  step. 
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Hamilton's  Principle 

The  method  used  in  this  section  for  deriving  first-order 
difference  equations  by  direct  application  of  Hamilton's  principle  was 
developed  by  Vance  and  Sitchin  [18].   The  section  on  the  derivation  of 
the  method  from  [18]  is  included  in  the  Appendix  for  completeness. 

The  motion  of  a  dynamical  system  is  determined  by  solving  the 
3nN  equations  derived  from 

oS     2J,      df..      ^ 

sr=-  +     L     X.  .  -r-^2-  -  "  F.,  At  ,  (4-10) 

^   J-l   1J  °?ik       lk 

as    N    5f .  . 

— —  +  £  X .  .  — =3.  =  o  ,  (4-ii) 

ap. _      j=i    1J  aj>. . 

ik  Tik 


f.,  =  Ap.,  -cp.,  At  =  o,  (i=l,2,...,n;  k=  1,2,...  ,N) ,    (4-12) 

ik     ik    ik 


where 


and 


\        -  i    undetermined  Lagrangian  multiplier  at  j    time 
ij 

interval , 

*     th 
F   -  i    generalized  equivalent  nonpotential  force,  at 

,  th   *  • 

k        time   interval, 

N 

S     =    Y     L*(c?      ,   cp      cp      ,    cp  cp      , .  .  .cp      )  At  .  (4-13) 

1        M       k^Ylk'   Y2k vnk'   Ylk'   Y2k'        Ynk' 

k=l 


The  3nN  first-order  difference  equations  resulting  from  Equa- 
tions (4-10)  -  (4-12)  in  conjunction  with  an  appropriate  difference 
expression  for  Up.      are  obtained  from  only  one  differentiation  and  are 
in  a  form  which  allows  their  solution  to  be  marched  out  with  time. 

Furthermore,  the  differentiations  (3S  /cf9 .)    and  (dS  /<&>   )  produce  the 

11         1    l 

.P   terms  and  uncouple  the  equations  in  one  step. 
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Example 

In  order  to  compare  the  two  methods,  the  above  procedure  will 
be  applied  to  the  system  in  Figure  4-1.  Since  the  internal  spring  K„ 
produces  an  internal  torque  which  is  transf errable  to  each  input  as  a 

generalized  equivalent  torque  (see  Chapter  II) ,  it  will  be  included  in 

* 
the  F   terms,  leaving  the  sum  S  made  up  of  only  the  kinetic  energy  of 

the  system.   Therefore, 

Sl  =  2  1  Zll  *1  +  2I12  V2  +  X22  *2  I  •  (4"14) 

Plk  ="(C11  91  +  °1*2  O  +  Ts/l,k-l  +  Vlk  *?lk  +  Tl'     (4"15) 

K 

F2k  =  "  (Si  *1  +  C2*2  O  +  TI/2,k-l  +  Kb/2k  *2k  +  T2'    (4"16) 

k 

A?   =  <p    -  9.,    ,  (4-17) 

xk    lk    ik-1 

where  the  C     and  K  ,   are  defined  in  Table  3-1. 
ij ,k      b/ik 

Substitution  of  Equation  (4-17)  into  Equation  (4-12)  reduces 
Equations  (4-10)  and  (4-11)  to  the  following  forms 

,9S  s 

/     1\  * 


and 


W  <<srl+  (Vk  •  (4-18) 

i   k 


»i.k"tr).  (4"19> 


Sep .   k 


Let   X.     be  given  by  the  finite  difference  relation 

X.  .  ,  =  (X.  .  ,  -  X.  .  )/At  .  (4-20) 

i,k+l     i,k+l    i,k 
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Then  Equation  (4-18)  becomes 

r,ds 


K.    ,  ,  =  A .  , 
1 ,  k+ 1    l ,  k 


w?.;k      i  kj 


(4-21) 


Substituting  Equations  (4-14)  -  (4-17)  into  Equations  (4-12) ,  (4-19) , 
and  (4-21),  and  solving  for  cp.     yields  the  following  set  of  six 


equations 


f 

K 

K 

l 

i 

= 

+  { 

1 
2 

_\ 

_\ 

k+l 

k 

i 

p      p      p      p 
11  1    12  1    21  1    22  1 


P       P      P       P 
11  2    12  2    21  2    22  2 


cpi  0 

92  0 

0  cp] 

0  cpr 


u^k 


-!k 


I-  ->k 


11 


'21 


* 
12 

»1 

+ 

T 
s/1 

+ 

* 
22_ 

k 

j2_ 

k 

T 
_  s/2_ 

k-1 

_   \ 


Kb/1  ^l 


Kb/2  ^2 


)  At,  (4-22) 


-Ikj 


9. 

1_  '-.to 

and 


L-  -Ik 


At 


(4-23) 


-1 


L-  Jk+1   L 


*11    X12 


X21    I22 


(4-24) 


-'k+l  "-  -'k+l 
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The  initial  values  of  the  X' s  are  obtained  by  solving  Equation  (4-19), 
using  the  initial  values  for  the  cp's  and  cp's.   Hence 


L>J0 


hi        '12 


X21    X22 


9. 


(4-25) 


The  matrix  Equations  (4-22)  -  (4-24)  require  one  less  numerical  step  of 
integration  as  pointed  out  in  I"30]-   It  is  also  pointed  out  in  refer- 
ence [30]  that  the  X' s  are  the  momenta,  usually  called  p.'s,  and  that 
Hamilton' s  canonical  equations  are  obtained  directly  without  formerly 
deriving  the  generalized  momenta.   These  equations  may  now  be  marched 
out  with  time  to   obtain  the  solution  to  the  dynamical  equations  of 
motion.   They  require  one  less  integration  step  as  opposed  to  solving 
Equations  (4-8)  and  (4-9)  numerically.   The  matrix  inversion  is  only 
necessary  once  per  integration  and  hence  any  error  is  only  integrated 
once  as  opposed  to  twice  for  Lagrange's  equations. 

The  significance  of  this  development  is  the  system  formulation 
in  terms  of  kinematic  influence  coefficients  developed  in  Chapters  II 
and  III,  allowing  complex  mechanisms  with  multiple  inputs  to  be  analyzed 
for  their  dynamic  response  by  established  numerical  methods.   This   for- 
mulation provides  a  way  to  reduce  complex,  multiple  input  mechanisms  to 
coupled  equivalent  mass  systems,  yielding  differential  equations  of 
motion  possessing  variable  coefficients.   These  coefficients  are  known 
in  terms  of  the  mechanism  geometry  through  the  kinematic  influence  coef- 
ficients of  velocity  and  acceleration. 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 


Complex  multiple  input  linkage  systems  have  been  difficult  to 
analyze  for  their  dynamic  response  because  of  their  nonlinear  geometric 
character.   This  characteristic  generates  holonomic  constraints  asso- 
ciated with  the  generalized  coordinates  necessary  in  describing  the 
motion  of  the  linkage  system.   Large  numbers  of  generalized  coordinates 
(see  Figure  5-1  (a))  have  been  required  to  obtain  the  system's  dynamical 
equations  of  motion.   Algebraic  equations  of  constraint  are  required  to 
account  for  the  generalized  coordinates  other  than  those  associated  with 
independent  system  inputs.   The  result  is  a  large  number  of  coupled, 
nonlinear,  second-order  differential  equations  together  with  a  set  of 
algebraic  equations  in  terms  of  undetermined  Lagrangian  multipliers 
which  account  for  the  geometrical  constraints  on  the  system.   The  algo- 
rithm required  to  integrate  numerically  and  solve  this  set  of  equations 
requires  a  dual  iteration  scheme,  one  to  solve  the  differential  equa- 
tions of  motion  and  one  to  satisfy  the  geometric  constraints  on  the 
linkage  system. 

Methods  such  as  those  developed  by  Chace  [11],  Uicker  [13], 
and  Carson  and  Trummel  [14],  utilizing  relative  coordinates  and  4x4 
matrix  coordinate  transformations,  have  been  the  only  tools  available 
to  solve  the  dynamic  response  question  for  these  systems.   The  set  of 
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N  =  8 

n  =  2 

8    Generalized    Coordinates 

6    Algebraic   Equations   of    Constraint 
(a) 


Figure  5-1.   Complex  Multiple  Input  Mechanism  and 
Its  Coupled  Equivalent  Mass  System 
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second-order  differential  equations  resulting  from  their  methods  has 
been  sufficient,  though  unwieldy,  time-consuming,  and  subject  to  error, 
to  describe  the  dynamical  behavior  of  linkage  systems. 

The  goal  of  this  dissertation  has  been  threefold:   to  develop  a 
systematic  method  whereby  linkage  systems  of  high-order-complexity  can 
be  constructed  from  Assur  groups  in  terms  of  kinematic  influence  coef- 
ficients of  velocity  and  acceleration  of  the  basic  system  group;  to 
reduce  these  highly  complex  linkage  systems  to  coupled,  equivalent  mass 
systems  acted  upon  by  equivalent  variable  rate  springs,  variable  coeffi- 
cient viscous  dampers,  and  equivalent  external  forces  and  torques;  and, 
to  determine  the  differential  equations  of  motion  for  the  coupled  equiv- 
alent mass  system  in  terms  of  the  minimum  number  of  generalized  coordi- 
nates (i.e. ,  the  number  of  independent  system  inputs). 

The  construction  of  general  linkage  systems  of  higher-order- 
complexity,  as  discussed  in  Chapter  II,  is  seen  to  be  expressed  in  terms 
of  series  and  parallel  link  connections.   The  connection  types  are 
defined  by  multiplication  (scries)  and  addition  (parallel)  of  successive 
velocity  influence  coefficients.   The  use  of  Assur  groups  to  construct 
mechanisms  of  higher-order-complexity  from  existing  system  groups,  with- 
out modifying  the  mobility  of  the  basic  chains,  allows  the  displacement 
analysis  to  be  performed  by  established  procedures.   This  is  the  basis 
for  eliminating  the  holonomic  constraints  on  the  system,  reducing  the 
number  of  generalized  coordinates  required  to  describe  the  motion  from 
N  (number  of  system  links)  to  n  (number  of  independent  system  inputs). 

Elimination  of  the  holonomic  constraints  subsequently  reduces 
complex,  multiple  input  linkage  systems  (see  Figure  5-l(a))  to  coupled, 
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equivalent  mass  systems  (Figure  5-l(b)),  requiring  only  n  generalized 
coordinates  to  describe  its  motion.   This  reduction  eliminates  the 
need  for  relative  coordinates  and  their  matrix  transformations  required 
by  the  existing  methods.   Rather,  the  coefficients  of  the  resulting 
differential  equations  of  motion  become  known  variables  of  the  system's 
independent  input  parameters  expressed  in  terms  of  kinematic  influence 
coefficients. 

Second-order  differential  equations  describing  the  dynamical 
behavior  of  the  equivalent  mass  system  have  been  derived  by  the  clas- 
sical Lagrangian  method,  while  first-order  difference  equations  were 
derived  by  the  direct  application  of  Hamilton's  principle.   This  method 
[18],  yields  a  set  of  first-order  difference  equations  derived  through 
only  one  differentiation.   These  difference  equations  can  be  marched  out 
with  time,  requiring  only  one  matrix  inversion  per  integration  step. 

The  equivalent  mass  system  formulation  developed  here  provides 
a  convenient  and  unique  medium  through  which  many  problems  concerning 
the  dynamical  behavior  of  linkage  systems  can  be  studied  and  simulated. 

The  influence  of  elastic  deformation  of  system  links  on  the 
dynamic  response  of  the  linkage  system  primary  input  can  be  studied, 
based  on  a  hinged  beam  model  such  as  the  one  shown  in  Figure  5-2. 
The  effect  of  bearing  deformation  on  the  input  dynamic  response  could 
be  investigated,  based  on  the  system  models  in  Figure  5-3.   The  equiv- 
alent mass  systems  are  shown  below  the  system  models. 

The  elimination  of  the  holonomic  constraints  places  complex 
linkage  systems  in  a  convenient  form  for  directly  applying  the  prin- 
ciples of  optimal  control  theory.   Optimal  open  loop  control  laws  may 
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Figure  5-2.   Linkage  System  with  Elastic  Coupler  Link 
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be  determined  which  minimize  specified  performance  indexes  while  satis- 
fying prescribed  constraints  on  the  control  itself,  the  states,  or  both. 
Figure  5-4  illustrates  this  concept  as  applied  to  complex  linkage  systems. 
The  problem  could  be  formulated  as  follows:   Determine  the  control  u 
which  minimizes  the  variation  of  the  velocity  s  from  the  velocity  s 

associated  with  position  s  over  the  range  s  £  S  £  S„,  subject  to  the 

1  \  & 

.* 

inequality  constraint  on  the  control,  u  ^  u  .   This  in  essence  places  an 

upper  bound  on  the  jerk  u,thus  giving  third-order  control  for  the  cam 

surface  designed  to  produce  the  required  control  u. 

The  preceding  problems  are  not  intended  to  be  solved  here. 

Rather,  they  are  provided  to  point  out  possible  research  areas  which 

can  be  pursued,  utilizing  the  system  formulation  developed  in  this 

dissertation. 
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APPENDIX 


APPENDIX  A 


DIRECT  DERIVATION  OF  FIRST-ORDER  DIFFERENCE 
EQUATIONS  FOR  DYNAMICAL  SYSTEMS 


The  derivation  of  the  set  of  first-order  difference  equations 
used  in  Chapter  IV  to  solve  the  time  response  of  complex  linkage  systems 
is  presented  in  part  from  the  paper  [18]     "Derivation  of  First-Order 
Difference  Equations  for  Dynamical  Systems  by  Direct  Application  of 
Hamilton's  Principle"  by  Vance  and  Sitchin.   The  purpose  for  presenting 
this  derivation  is  for  completeness  and  convenience  to  the  reader  of 
this  dissertation,  since  the  method  proves  to  be  ideal  for  treating 
complex  multiple  input  linkage  systems.   The  nomenclature  in  this  der- 
ivation does  not  correspond  directly  with  that  in  the  main  text.   It  is 
therefore  listed  separately  at  the  end  of  this  appendix. 

Derivation  of  Method 

Hamilton's  principle  for  nonconservative  systems  with  k 
degrees  of  freedom  is 

T       T   ^ 

r  Ldt  +  j* 

0        0   '1=1 


T        T  /  »        . 

6l  =  6  J*  Ldt  +  J*   (£  F.Sqldt  =  0,  (A-l) 


where  the  integrand  of  the  second  integral  is  the  virtual 
work  of  the  nonconservative  forces.   This  second  integral  is 
zero  for  conservative  systems.   After  partitioning  the 
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interval  0  to  T  into  N  small  increments  AT,  the  two  integrals 
can  be  approximated  by  sums  and  the  principle  can  be  rewritten 
as 

6S1  +  S2  =  0,  (A-2) 

where  the  functions  S  and  S  are  sums  given  by 
N 

Sl  =  E  Ln(qin'q2n'---'qkn'V'V'"-'\n)At  '     (A"3) 
n=l 

N    k 

S2  -  E   E  Fin6qinAt'  (A"4) 

n=l  i=l 

Equation  (A-2)  requires  that  the  variation  of  the  function 
S  equal  the  negative  S  .   Since  the  displacements  and  veloc- 
ities are  to  be  related  by  some  finite-difference  expression, 
they  are  not  independent  of  each  other.   The  problem  lends 
itself  to  the  use  of  Lagrangian  multipliers  in  order  to  achieve 
an  independent  variation  of  coordinates.   The  equations  of 
constraint,  defined  by  the  previously  mentioned  finite- 
difference  relationship,  have  the  general  form 

Aq. 

in  . .  _« 

in    At 


f.   =  Aq.   -  q.  At  =  0,  (A-S) 

in     in    m 

where  Aq.   is  any  desired  expression  for  the  first-order  dif- 
ference of  q.  . 
in 
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Taking  the  variation  of  the  constraint  function  (A-6) 
gives 


N   /0f 


6f-   =  E 

in     *—> 


df. 
6q.  .  +  -rill   6q.  .  )  =  0. 


3=1  V^ij  "iJ  '  6^ij   ^ 


(A-7) 


Equation  (A-2)  can  now  be  written  as 

N   k 


6S 


,  +  I!    Z  *.  6f-   +  s„  =  o. 

1    *-*,  .V,   in  in    2 


n=l  i=l 


(A-8) 


Substitution  of  Equation  (A-7)  into  Equation  (A-8)  and 
rearrangement  of  terms  gives 


N 

k 

E 

E 

n=l 

i=l 

I 

n  df. 


-     +  F.  At  +  V     \.  .   5~i 

in      .«   ij  oq. 
in  j=l       in 


6q. 
in 


N   k 

E  E 

n=l  i=l 


3*7  +  A  *» 


of.  . 
ij 


in  j=l   1J  °qin 


6q.   =  0. 
in 


(A-9) 


The  kN  X   are  chosen  so  that  the  kN  bracketed  expressions 
ij 

in  the  second  double  summation  are  all  zero.   This  leaves  the 

6q    as  independent  variations.   Equation  (A-9)  can  then  be 
in 

satisfied  by  independently  requiring  the  kN  bracketed  expres- 
sions in  the  first  double  summation  to  be  zero.   Thus  the 
motion  of  the  dynamical  system  will  be  such  that  the  following 
3kN  equations  are  satisfied  (i  =  l,2,...,k;  n=  1,2,...,N): 
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as,  af .  . 

^   +  ^  Vj  ^   =  -  F.  At,  (A- 10) 

in   j  =  l       in      in 

3S.     N      Bf .  . 

557  +  ?,  x«  **  ■  "■  (A"11} 

in   j=l       in 


f .   =  Aq.   -  q.  At  =  0.  (A-12) 

in     in    in 


When  the  first-order  difference  form  Aq.   is  substituted 

in 

into  Equation  (A-12) ,  the  summations  in  Equations  (A-10)  and 
(A-ll)  are  reduced  to  only  a  few  terms.   For  example,  if 

Aq .   =  q .   -  q .  '■  ,  , 
in    in    in-1 

3f .  .  of .  . 

3-22  _  o,  j  4   n,  n  +  1;  and  -^^-  =0,  if  j  4   n.    (A-13) 

in  in 

As  is  characteristic  of  the  Lagrangian  multiplier  method, 
the  convenience  of  treating  dependent  variables  as  if  they 
were  independent  has  been  gained  at  the  expense  of  an  added 
set  of  unknowns,  the  \'s.   Unlike  many  applications  of 
Lagrangian  multipliers,  however,  the  \' s  are  not  in  general 

constant.    In  fact  the  X' s  represent  momenta,  .  .  .,  and  have 

** 
the  status  of  independent  coordinates  [3]  . 

The  set  of  Equations  (A-12)  may  be  considered  trivial 
(although  necessary)  in  the  sense  that  they  are  simply  the 
equations  of  constraint  between  the  velocities  and  displace- 
ments. 


The  \' s  will  be  constant  only  in  the  case  of  ignorable 
coordinates. 

Goldstein,  H. ,  Classical  Mechanics,  Addison-Wesley  Publishing 
Co.,  Atlanta,  1965,  p.  227. 
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The  derivation  just  presented  is  not  tied  to  any  particular 
finite-difference  form.   A  determination  of  the  best  form  to 
use  will  depend  on  the  application. 

The  symbols  are  defined  as 

6  =  "variation  of" 

A  =  finite  increment 

.  =  dot  appearing  directly  above  a  variable  designates 

derivative  with  respect  to  time 
X  =  undetermined  Lagrangian  multiplier  constant 
i  =  the  corresponding  generalized  coordinate 
n  =  the  corresponding  time  interval 
j  =  the  corresponding  time  interval 
t  =  time 
F.   =  generalized  nonpotential  forces 


APPENDIX  B 


NUMERICAL  SOLUTION  TO  A  TWO  DEGREES 
OF  FREEDOM  EXAMPLE 


The  purpose  of  this  appendix  is  to  illustrate  the  actual 
responses  obtained  from  the  dynamical  equations  of  motion  derived 
for  the  five-bar  mechanism  in  Chapter  IV  (see  Figure  (B-l(a))). 

The  procedure  for  obtaining  the  variable  coefficients  to 
Equations  (4-22)  -  (4-25) ,  describing  the  dynamical  behavior  of  the 

system  in  Figure  B-l(a),  is  described  below.   The  necessary  coefficients 

*      *      *     *  * 

are  K,  ...  ,  C.  .  .  ,  T.  ,  ,  I.  .  .  ,  and  .  .P  ,  . 

X/ik   xj,k   i,k'   ij,k'     ij  r,k 

(1)   For  the  given  initial  position  of  the  mechanism  defined 


by  cp  and  cp   calculate  the  kinematic  position  informa- 
tion pertaining  to  the  other  links  §  ,  $  ,  etc. ,  by 
using  the  equations  in  Figure  B-2. 

(2)  Substitute  this  position  information  into  the  appropriate 
blocks  of  Table  2-1  to  determine  the  velocity  and  acceler- 
ation influence  coefficients  ..g.  and  „h...   For  example, 

Jfc  1      X  ij  r       , 

$      _  $ 
3  l,k+l   3  l,k-l 
3gi  -  25jT ' '  CB-1) 

(3)  Once  the  »g.'s  and  „h..'s  are  determined,  substitute  them 

*  l        x  IJ 

into  the  appropriate  blocks  of  Tables  3-1  and  3-2.  This 
furnishes  the  following  expressions  for  the  coefficients 
to  the  differential  equations: 
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2         2 

*        N   /_  \ 

c  •  „  =  I  (c,(,g.)(,e.)J  .  (b-3) 

1 J  '    jtl       1    j  k 

T*    =  (,B.„)T,  •  (B-4) 

xk      Z  lk   £ 

I     and   P    are  determined  from  the  expressions  in  Table  3-2, 
ij,k     ij  r,k 

blocks  1.1  and  1.2.   The  right-hand  side  of  Equations  (4-22)  -  (4-25)  are 

st 

now  completely  known,  allowing  the  k+1   \'s,  cp's  and  cp's,  to  be  deter- 
mined.  This  procedure  is  repeated  until  the  equations  of  motion  are 
integrated  over  a  predetermined  time  interval.   Equations  (4-22)  - 
(4-25)  represent  the  equations  of  motion  for  the  coupled  equivalent 
mass  system  shown  in  Figure  (B-l(b)). 

The  solutions  to  the  dynamical  equations  of  motion  are  shown 
in  Figure  (B-3)  and  (B-4).  Polar  plots  of  <p  vs  cp  and  <p  vs  cp  are 
shown  for  the  parameters  listed  in  Table  B-l. 

For  added  clarity,  Figure  (B-5)  illustrates  the  equivalent 
inertias  as  functions  of  time.   The  equivalent  inertias  show  no  cyclic 
phenomenon  due  to  the  noncyclic  character  of  the  input  links.   However, 

the  inertia  coupling  term,  I   ,  illustrates  the  coupling  between  the 

*       * 
positive  inertia  terms  I   and  I22- 

The  influence  coefficients  were  calculated  by  finite  differences. 
The  total  problem  was  programmed  on  the  IBM  360-65.   The  program  con- 
sisted of  270  cards  and  required  0.38  minutes  to  execute  for  an  integra- 
tion step  size  of  At  =  .000025  second.   This  step  size  can  be  increased 
considerably  without  affecting  the  solution  accuracy,  thus  decreasing 
the  execution  time  by  that  factor. 
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(a) 
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(b) 


Figure  B-l.   Two  Degrees  of  Freedom  Five-Bar  Example 
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/34  =  A  +  A    -    2/3Acosc£ 


34 


/34=  (./5COS(£6  -/|C0Sc/>|+i'2C0S<^2+(/5SIN(^5+  ^SIN^>r/2SIN^ 


*34=      C0S 


-i 


•*3      -*4     ^34 


$-  s      SIN 


-I 


2AA 

^5SlN(/>5+   ^SINcjb,-    /2SIN<^>2 


34 


$4*     SIN* 


("  ^5SIN^>3+  -^,  SIN</> ,  "    ^SlNC^j 


34 


+  SIN 


SIN 


Figure  B-2.      Kinematic  Position  Equations 
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TABLE  B-l 
FIVE-BAR  PARAMETERS 
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Length  (inches) 

Weight  (lbs. ) 

Moment  of  inertia 
(in.  lb.  sec.2) 


4.0 

2.0 

8.0 

12.0 

11.18 

2.97 

1.48 

5.52 

4.16 

-  0  - 

.0488 

1.14 

.1385 

.4062 

-  0  - 

K„,«  =  50.0  in.  lb./rad. 
34 

C   =   5.0  in.  lb.  sec./rad. 
T   =    .05  in.  lb. 
T2  =   -.02  in.  lb. 


(«Vo 


Op..)    =  250  rad./sec. 


(Cf>  )    =  90c 
Z    o 


(Vo   "   ° 


Solution  to  differential  equations  for  known  forcing 
functions,  constant  T.  and  T  ,  are  shown  on  Figures 
(B-3)  and  (B-4).     l  l 
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